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A new quasilinear mathematical model of heat conduction with finite velocity 
of heat front movement is offered. 

i_i , 1. The analysis of existing models of heat conduction. 

^5 \ Now as basic mathematical model in the theory of heat conduction the linear 

lO ' parabolic differential equation is used. In case of one spatial dimensional: 
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t^j- ' Also it can be presented as system of two equations: 
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O 1 °P^7 = -dwJ, (2) 
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j = -a (3 ) 

^ The following designations are entered: T(t, x) - field of temperature, J(t, x) - 

vector of heat flow, A - coefficient of heat conduction, Dt = A/(cp) - coefficient 
of temperature conduction, p - density, c - thermal capacity, div = d/dx, i - unit 
vector. 

, ^ ' Equation (1) is obtained taking into account suppositions about instantaneous 

^ ■ relaxation of the heat flow, that contradicts physical conceptions about the nature 

of heat. This leads to infinite velocity of front heat movement in the given model. 
Sometimes model the nonlinear heat conduction equation: 

dT d dT 

where the coefficient of heat conduction is function of temperature, or linear 
hyperbolic equation Cattaneo is used. But both these models have basic lacks[2], 
which do not allow to widely use them at the decision of practical tasks. 

In this work new heat conduction equation is obtained, which is a generaliza- 
tion of the linear equation (1). For this purpose equation (1) is transformed to 
the nonlinear equation in such a manner, that quantitatively its solutions differ 
from solutions of equation (1) a little, but the new equation will describe the 
process of heat conductivity qualitatively precisely. 

In the new equation there is a parameter, connected with velocity of heat 
conduction front. 
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2. The new model of heat conduction and investigation of its local 
groups of symmetry. 

Let's change the equation (3) and have it presented as 

J = - l X{—-a.arctg{-—)), (5) 

where a - parameter, which has dimension of a gradient of temperature. From 
equation (5) it is clear, that flow J on the absolute value never exceeds the flow 
described in equation (3) more than 0.5a7r and relative value of the flows differ- 
ences approaches zero, when temperature gradient increases tending to infinity 
or this gradient being close to zero. Equation (5) and equation (2) will give the 
new equation of heat conduction: 

or {dr/dxf o 2 t 

dt T {dT/dxf + a 2 dx 2 ' [ ) 

Equation (6) is invariant concerning the local groups of symmetry: 

x 1 = d t , x 2 = d x ,x 3 = d T , x A = 2td t + xd x + Td T . (7) 

Using invariants T/y/i = const and xj\Jt = const generated by the operator 
of dilatation X^ it is possible to obtain replacement of variables z = x/\/t, 
T(x,t) = \/if(z). Substituting them into equation (6), we obtain the ordinary 
differential equation for function f(z) : 

3. Numeric solution of equation (8). 

Consider boundary value problem for equation (6) with initial condition T(t = 
0,x) = 0, boundary conditions T(t,x = 0) = By/i,dT(t,x = 0)/dx = C, x > 
and assume, that Dt = 1- 

At reduction of equation (6) to equation (8) boundary conditions will be 
transformed to a kind of /(0) = B,df(0)/dz = C. In figure 1 results of the 
numerical solution of the equation (8) at B = 1, C = —1, a 2 = 0.001, z = 
0...5 are shown. Velocity of movement of heat front being Vq = Vo(i) it is possible 
to obtain the following: 

f(z = Zq) = — >■ z = xo/Vt — >■ Vq = dxo/dt = z / (2\/i). 

Function / = f(z) accepts zero value at z ~ 1.55 and, consequently, velocity of 
movement of heat front will be Vq ~ 1.55/ (2\/t). 

This solution qualitatively precisely describes the process of heat conduction 
in real physical media. 

4. Example of non-local symmetry of the equation (6). 

Equation (6) by means of non-local transformation H(t,x) = dT/dx can be 
transformed to a kind of: 

dH__ d_ ( D T {dH/dx) 2 dH \ 
~dt~dx \ (dH/dx) 2 + a 2 dx ) ' 
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This is a nonlinear equation of heat conduction, in which the coefficient of heat 
conduction is function from gradient of temperature. 
5. Quasilinearity of equation (6). 

Numerical researches of equation (8) have shown, that its solutions have phys- 
ical sense at < a < 0.1. It is possible to assume, that this parameter and in 
the initial equation (6) also will accept only such small values. Then equation 
(6) can be presented as a system of the equations: 

dr n d 2 T c (d 2 r\ 2 2 

1h™ Dt 1^> for [d?) >>a ' (9) 

dT d 2 T 

-.0, for^O. (10) 

Equations (9) and (10) are close to a linear kind and, therefore, equation (6) can 
be named as the quasilinear equation of heat conduction. 

l.Tikhonov A.N., Samarsky A. A. Equations of mathematical physics. - Moscow: Nauka, 1966. 
- 660p. (in Russian) 

2. Lykov A. V. Theory of heat conduction. - Moscow: Vyshaya shkola, 1967. - 559p. (in Russian) 



«*5i 




3 



